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Abstract 

The scaling ranges of time correlations in the cloud base height records of 
marine boundary layer stratocumulus are studied applying the Detrended 
Fluctuation Analysis statistical method. We have found that time dependent 
variations in the evolution of the a exponent reflect the diurnal dynamics 
of cloud base height fluctuations in the marine boundary layer. In general, a 
more stable structure of the boundary layer corresponds to a lower value of the 
a - indicator, i.e. larger anti-persistence, thus a set of fluctuations tending to 
induce a greater stability of the stratocumulus. In contrast, during periods of 
higher instability in the marine boundary, less anti-persistent (more persistent 
like) behavior of the system drags it out of equilibrium, corresponding to larger 
a values. From an analysis of the frequency spectrum, the stratocumulus 
base height evolution is found to be a non-stationary process with stationary 
increments. The occurrence of these statistics in cloud base height fluctuations 
suggests the usefulness of similar studies for the radiation transfer dynamics 
modeling. 
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I. INTRODUCTION 



Cloud base height (CBH) profiles are difficult object to describe and analyze. The base 
of a cloud is determined by the local condensation level up to which unsaturated air parcels 
rise. CBH has a complex structure, that varies with the coordinates of the observation 
point and time. Such a complexity is due to a variety of processes that take place in the 
marine boundary layer (MBL) formed over the sea surface. Some of these processes are, e.g., 
turbulent motions, entrainment,Bi radiative transfer and changes in the cloud microphysical 
structure. Thus the CBH looks like an erratic quantity in time and space. Therefore, 
advanced methods of computational statistical analysis are required in order to extract 
meaningful physical information from data series with such a high level of complexity that 
traditional analysis techniques fail to provide. The detrended ffuctuation analysis (DFA) 
method is designed to search for long-range correlation^ in non- stationary, highly ffuctuating 
signals. The method has been successfully applied in the investigation of turbulence, 
sequences,^ ionic transport trough cell membrane,Q'Ll foreign currency exchange rates, 
meteorological phenomena,inhEi and stratus clouds dynamics.0'El Using this techniquelli 
we hereby sort out correlations and anti-correlations in stratocumulus cloud base height 
ffuctuation data records. Our ffndings of the time evolution of the local a-exponent (which 
is the Hurst exponent^), i.e. describing a scaling law over a certain time range, are related 
to the physical processes in the cloud. These results reffect the ability of the method to 
indicate the type of physics that underlies the CBH profile phenomenon. 



II. MARINE BOUNDARY LAYER - STRUCTURE AND PROCESSES 

A variety of physical processes takes place in the atmospheric boundary layer. At time 
scales of less than one day, significant fiuxes of heat, water vapor and momentum are ex- 
changed due to entrainment, radiative transfer, and/or turbulence.lli@ The turbulent char- 
acter of the motion in the atmospheric boundary layer (ABL) is one of its most important 
features. The turbulence^ can be caused by a variety of processes, like thermal convections, 
or mechanicallvgenerated by wind shear, or following interactions infiuenced by the rotation 
of the Earth.EiE3 This complexity of physical processes and interactions between them cre- 
ate a variety of atmospheric formations. In particular, in a cloudy ABL the radiative fiuxes 
produce local sources of heating or cooling within the mixed-layer and therefore can greatly 
infiuence its turbulent structure and dynamics. Moreover the variations in the turbulent 
structure and dynamics of the clouds cause subsequent changes in the cloud boundaries, 
especially in the cloud base. 

The atmospheric boundary layer is characterized by an inner (surface) layeiSJ^l at heights 
above the aerodynamic roughness length@ and below one tenth of the depth of the ABL, 
while the outer (Eckman) layer is at higher distances within the ABL. The outer region of an 
unstably stratified ABL is often referred to as the mixed layerSS because of the dominating 
convective motions that take place there, generated by strong surface heating from the Sun 
or by cloud-topped radiative cooling processes.0 In contrast, a stably stratified ABL occurs 
mostly at night in response to the surface cooling due to long-wave length radiation emitted 
into the space. 
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The top of the boundary layer in convective conditions is often well defined by the 
existence of a stable layer (capping or subsidence inversion), beneath which clouds form, the 
so-called cloud-topped boundary layer (CTBL).il In presence of clouds (shallow cumulus, 
stratocumulus (Sc) or stratus (St)) the structure of the ABL is modified because of the 
radiative fluxes. Phase changes become more important. During cloudy conditions one can 
distinguish mainly : (i) the case in which the cloud and the subcloud layers are fully coupled; 
(ii) two or more cloud layers beneath the inversion, with the lower layer well-mixed with 
an upper elevated layer, decoupled from the surface mixed layer or (iii) a radiatively driven 
elevated mixed cloud layer, decoupled from the surface. 

An interesting case is that of the marine boundary layer characterized by high concentra- 
tion of moisture. It is wet, mobile and has a well expressed lower boundary. The competition 
between the processes of radiative cooling, entrainment of warm and dry air from above the 
cloud and turbulent buoyancy fluxes determine the state of equilibrium of the cloud-topped 
marine boundary layer .t3 The Atmospheric Stratocumulus Transition Experiment (ASTEX) 
was designed to clarify the transition from stratocumulus to trade cumulus cloudsEl in the 
MBL in the region of the Azores Islands.S Several papers have been dedicated to the rel- 
evance of the main processes causing this transition.Sil In Ref. [27] Betts et al. found a 
relationship between the thermodynamic structure of the marine boundary layer and the 
diurnal variation of the cloudiness. 



III. DFA METHOD AND SPECTRAL ANALYSIS 

A method that relaxes the requirement of stationarity of the investigated signal is the 
detrended fluctuation analysis (DFA) method. The DFA method is a tool used for sorting 
out long range correlations in a non-stationary self-affine time series with stationary incre- 
ments. The method has been used previously in the meteorological field.00'0 It provides 
a simple quantitative parameter - the scaling exponent a, which is a signature of the cor- 
relation properties of the signal. The Detrended Fluctuation Analysis technique consists^ 
in dividing a random variable sequence y{n) of length into N/t (non-overlapping) boxes, 
each containing t points {N/t = 4,5, ...). Then, the trend (assumed to be linear in this 
investigation, but it can be generalized without any difficultyB) z(n) = an + b in each box 
is computed using a linear least-square fit to the data points in that box. The detrended 
fluctuation function F'^{t) is then calculated following : 

1 / /V \ 

F\t) = - J: [yin) - z{n)]' for A; = 0, 1, 2, . . . , - - 1 . (1) 

Averaging -F^(t) over the N/t intervals gives the mean-square fluctuations 

< F^t) >i/2_ (2) 

If the y{n) data are random uncorrelated variables or short range correlated variables, the 
behavior is expected to be a power lawS with an exponent^ a=l/2 if the fluctuations are 
not correlated. An exponent a ^ 1/2 in a certain range of t values implies the existence of 
long-range correlations in that time interval as, for example, in fractional Brownian motion.E2l 
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A small value of a indicates antipersistenceEj of correlations, i.e. a positive fluctuation is 
more likely to be followed by a negative one than a positive one. An a value between 0.5 and 
1 indicates persistence of correlations, i.e. a positive fluctuation is more likely to be followed 
by a positive one than a negative one. The a-exponent value that holds true for a certain 
time interval called scaling range, is a characteristic of the correlations in the fluctuations of 
a signal y{t) defined between the beginning to and end of the observations tM, i-e. [tc^A/]- 
The value can vary in time if the process is not stationary and/or the scaling range is finite 
(see Sect. 6). The a exponent is related to the usual fractal dimension@'i3 D. 

The concept of the invariance of the fractal shapes^ at different scales (self-similarity) 
can be transposed in the statistical analysis of a time series. One searches for a eiven scale 
range where the statistical properties of the signal are invariant. By definitionE9 the time 
series is self-affine if its spectrum has a power law dependence. 

s{f) - {1/fr. (3) 

It has been shown by Heneghan and McDarbyil that the relationship P = 1 + 2a holds 
true for stochastic processes, i.e. for fractional Brownian walks. Depending on the value of 
the spectral exponent /3, one has the following caseJUHl : if /? < 1, the process is stationary; 
if /3 > 1, the process is non-stationary; if 1 < /? < 3 the process is non-stationary with 
stationary increments. In the case of white noise (3 = while (3 = 2 for a Brownian walk 
signal. 



IV. DATA 

Cloud base height data has been obtained during the Atlantic Stratocumulus Transition 
Experiment (ASTEX) in June 1992, at the Azores Islands using a laser ceilometer. The 
wavelength of the vertically transmitted into the atmosphere laser beam is 0.904 mm. The 
ceilometer works with a 15 m spatial resolution and reaches the maximum measurable height 
of 4 km. The time resolution of CBH records is 30 seconds. The cloud base height values 
were selected from the profiles of the backscattered signal intensity, using manufacturer's 
algorithm based on a voltage threshold for the maximum value of the backscattered radiation 
at a heightEj giving the local condensation level up to which unsaturated air parcels rise. 
We use data records measured on June 15, 1992 (Fig. 1(a)) and June 18, 1992 (Fig. 2(a)) 
to determine (i) the scaling properties of CBH fluctuations; (ii) to classify the type of time 
correlations, (iii) to investigate the diurnal changes in the CBH evolution, and (iv) to identify 
the type of the time dependent long-range correlation evolution. 

According to,0 on June 15, 1992 an anticyclone was centered north of Santa Maria - 
one of the Azores Islands, where the measurements were taking place, producing a strong 
subsidence inversion beneath which were observed marine boundary layer clouds. The noc- 
turnal MBL was found to be well-mixed, generally coupled to the surface, with a solid layer 
of stratocumulus, while the daytime MBL was observed to be decoupled from the surface, 
because of the coexistence of a stratocumulus cloud layer and a cumulus marine subcloud 
layer. After the sunrise the MBL became intermittently decoupled and CBH varied dra- 
matically as a cumulus formed at the top of the marine subcloud layer. In this decoupled 
MBL, the cumulus convection, that arises from the latent instability generated in the marine 
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subcloud layer, acts toward a recoupling of the formally decoupled cloud layers that persist 
beneath the trade inversion. Five sequences of strong decoupling were observed that day 
(Fig. 1(a)). 

For June 18, 1992 we do not have any information about the meteorological situation, 
but in analogy with data in Fig. 1(a), the CBH profile on Fig. 2(a) leads one to distinguish 
periods of increased variability of CBH, followed by short CBH fluctuations. 

V. SCALING PROPERTIES OF THE CBH DATA 

The CBH evolution, reported in Figs. 1(a) and 2(a) through S{f), can be regarded 
as a non- stationary process with stationary increments, since their spectral exponents are 
(3 = 1.28 ± 0.1, for June 15, 1992 (Fig. 3) and f3 = 1.49 ± 0.08, for June 18, 1992 (Fig. 4). 
The data error bars are calculated from the r.m.s deviations in the fit according to standard 
techniques and have their 95% confidence interval conventional meaning.ii 

The type of correlations of the CBH fluctuations can be probed using the Detrended 
Fluctuation Analysis (DFA) method. Applying the DFA method, we find that long-range 
time correlations in CBH fluctuations exist and are of an anti-persistent type with a = 
0.24 ± 0.002 for June 18, 1992 and a = 0.21 ± 0.005 for June 15, 1992 (Fig. 5). They 
hold for a time interval approximately equal to 140 min in the first case and 40 min for 
the second. In this scaling range, the relationship j3 = 1 + 2a is reasonably satisfied. The 
flat like regions of both curves (Fig. 5) correspond to low values of the Hurst exponent a 
(a = 0.06 ± 0.003 for June 15, and a = 0.10 ± 0.003 for June 18) and indicate the presence 
of 1/ /-like noise. 

For completeness, we show in Fig. 6, the distribution function of the CBH fluctuations for 
the 30 sec time lag interval. Notice the similarity between both cloud cases and the marked 
departure from the Gaussian behavior, as numerically exemplified in the S{f) and DFA 
correlation methods in fact. A similar double triangular pyramid shape has been observed 
in other meteorological cases.il 

VI. TIME CORRELATION ANALYSIS 

It is of interest to test if the correlations remain the same for shorter intervals accommo- 
dated into [^c^m] or if they change with time, as should be anticipated for nonstationary 
time series data. In order to probe the existence of so called locally correlated and decorre- 
lated sequences, one can construct an "observation box" with a certain width, u, of size 7h, 
place the box at the beginning of the data, calculate a for the data in that box, move the 
box by Auj= 30 min (or 60 points) toward the right along the signal sequence, calculate a 
in that box, a.s.o. up to the A^-th point of the available data. A time dependent a exponent 
may be expected, and is given at the end of each box, thus for t ranging from Au to N. 

Thus we obtain an "instantaneous measure" for the degree of long-range correlations. 
Results for this instantaneous a-exponents are shown in Figs. 1(b) and 2(b) for the cloud 
base height data on June 15, 1992 and June 18, 1992. Note that each a value represents the 
behavior resulting from all points in the box, which continuously overlap. This approach 
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seems suitable to cloud base height data because it can be expected to reveal changes in the 
correlations dynamics of the clouds at various times for a given time lag Au;. 

The a evolution is within the limits of a quite anti-persistent behavior with a transition 
from weak to strong, at various instant, in other words, anti-persistence is found in the 
underlying dynamics. Five regions in the CBH evolution, corresponding to the description 
given in section 2, can be distinguished^ for the CBH data measured on June 15. The first 
region lasts from to 6 h, the second from 6 to 13 h, the third holds from 13:00 to 16:00 h, 
the fourth - from 16:00 to 20:00 h and the fifth one from 20:00 to 24:00 h. The transition 
occurring from 08:00 to 10:30 in the evolution of the a exponent indicates the increasing 
role of a process that generates the second structure in the CBH that appears between 10-13 
h in Fig. 1(b). This process is characterized by a mean value of the local a = 0.24 ± 0.01 
and is understood as a process of destabilization, in the decoupled MBL. Between 14-18 h 
the downward trend of the a exponent indicates the increasing dominance of a process that 
generates the third part of CBH evolution that exists between 17-19 h in Fig. 1(b). The 
mean value of the local a exponent is a =0.20±0.02. From 19 h to 20:30 h a very slow 
upward trend leads to the process of the fourth region between 21:30 h - 24 h in Fig. 1(b), 
with a mean local a =0.26±0.02. 

For June 18 we distinguish three main regions in the CBH evolution: 0-10 h, 10-19 h and 
19-24 h. The trends in the a exponent evolution are as follows: 0-10 h - process, generating 
the first CBH structure, with mean local a =0.17±0.02. This corresponds to the dynamics 
of the MBL before its decoupling. For the period 10-16 h a large upward transition leads 
to a weak anti-persistent correlation in CBH fluctuations with a mean value of the local 
a =0.35±0.02. Hence the dominant process is destabilizing the MBL structure. For the 
time 16-20 h the dynamics of the process changes again, characterized by a mean local a 
=0.26±0.01 and a return to more stable MBL structure. 

Thus the evolution of the Hurst exponent a in the cases that we have investigated, con- 
sists in the successive transitions from lower to higher values and vice versa. This corresponds 
to transitions from strong to weak anti-persistent correlations and reflects the changes in 
the underlying dynamics of the MBL. 

Notice that the DFA function being a measure of the root mean square deviation from 
the linear best fit (a "local trend") is therefore a measure of a height-height correlation 
function. The so-called height is nothing else than a mass threshold. Therefore this DFA 
function informs on a corresponding physical property, i.e. the local compressibility (or bulk 
modulus) of the cloud fiuid.i3ii 

VII. CONCLUSIONS 

We have presented an analysis of a marine stratocumulus cloud base height data series 
using the DFA method and the spectral density method in order to investigate the scaling 
range of time correlations in the CBH fluctuations searching for scaling laws and persistence 
effects. In all cases anti-persistent correlations hold with a well established day-night evolu- 
tion measured by the instantaneous a index , - itself characterizing the dominant physical 
process. The transition, revealed by the trend in the a evolution, between processes in- 
ducing different levels of anti-persistent correlations of the fluctuations could be used as a 
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trace parameter of the CBH dynamics. The accuracy of these statistics in CBH fluctuations 
suggests the usefulness of similar studies for the radiation transfer dynamics. 

The DFA method allows one to obtain the scaling ranges and crossovers of the correlations 
between the fluctuations of a signal. The a exponent is in fact directly related to that for 
the auto correlation function.! Incidently, if the relationship 7 = 2(1 — a) is supposed to 
hold for describing the power law dependence of the autocorrelation function, a value 7 = 
1.56 is expected. This means that the CBH fluctuations can be empirically represented by 
a fractional Brownian motion process.!^ The oscillations between different stability regimes 
could thus be mapped onto a Kramer problem,!! itself mimicking the basic Langevin equation 
for the signal increments. 

Nevertheless this knowledge of the second moment of a signal is not usually enough in 
order to have a full understanding of the physical process nor of its underlying physical 
mechanisms governing the CBH. It would be of interest to couple these results to some 
analysis of the local (vertical and horizontal) velocity fluctuations, and of the temperature 
fluctuations in order to obtain additional information, whence the viscosity and thermal 
conductivity, for a better understanding of cloud physics and weather predictability. 
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FIG. 1. (a) CBH fluctuations for June 15, 1992; (b) The a evolution, clarifying the transitions 
between strong and weak anti-persistent behavior of CBH fluctuations. 
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FIG. 2. (a) GBH fluctuations for June 18, 1992; (b) The a evolution, clarifying the transitions 
between strong and weak anti-persistent behavior of GBH fluctuations. 
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FIG. 3. The power spectrum of the cloud base height data (data in Fig. 1(a)). A spectral 
exponent (3 =1.28±0.1 characterizes the correlations of fluctuations. 
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FIG. 4. The power spectrum of the cloud base height data (data in Fig. 2(a)). A spectral 
exponent /3=1.49±0.08 characterizes the correlations of fluctuations. 
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FIG. 5. The DFA function < F{t) >^/2 foj. CBH data measured on June 15, 1992 (data 
in Fig. 1(a)). Scaling properties change from ai = 0.21 ± 0.005 to 02 = 0.06 ± 0.003 for time 
lags longer than 40 min. CBH data measured on June 18, 1992 (data in Fig. 2(a)) scale with 
ai = 0.24 ± 0.002 for time lags smaller than 140 min and 0:2 = 0.10 ± 0.003 after that. 
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FIG. 6. Partial distribution function (PDF) of the CBH signal increments (discretization step 
= 30 s) for both June 15 (black dots) and June 18, 1992 (black squares) stratocumulus clouds 
observed on Azores Islands. The x-axis is normalized in units of the standard deviation of the 
corresponding PDF for a time lag 960 s, - supposedly long enough such that each PDF then 
corresponds to a Gaussian distribution. The standard deviation for June 15 is 198.9 m and the 
one for June 18 is 323.4 m. The June 15 data has been displaced vertically by one decade. 
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